Vortex core states in superfluid Fermi-Fermi mixtures with unequal masses 
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We analyze the vortex core states of two-species (mass imbalanced) superfluid fermion mixtures as 
a function of two-body binding energy in two dimensions. In particular we solve the Bogoliubov-de 
Gennes equations for a population balanced mixture of ^Li and ''"K atoms at zero temperature. 
We find that the vortex core is mostly occupied by the light-mass (^Li) fermions and that the core 
density of the heavy-mass (*°K) fermions is highly depleted. This is in contrast with the one-species 
(mass balanced) mixtures with balanced populations where an equal amount of density depletion is 
found at the vortex core for both pseudospin components. 

PACS numbers: 03.75.Hh, 03.75.Kk, 03.75.Ss 



Recent observation of quantized vortices in two- 
component mass and population balanced fermion mix- 
tures of ^Li atoms with short-ranged attractive interac- 
tions [l| has not only complemented the previously found 
evidence [1, H, 0, d, @| , but also indicated a very strong 
evidence for the superfluid ground state. In these exper- 
iments it has been realized that the ground state evolves 
smoothly from a Bardeen-Cooper-Schrieffer (BCS) su- 
perfluid to a molecular Bose-Einstein condensate (BEC) 
as the attractive interaction strength varies from small to 
large values, marking the first demonstration of the the- 
oretically predicted BCS-BEC crossover 0, 1, More 
recently quantized vortices have also been studied in 
mass balanced but population imbalanced fermion mix- 
tures [l0[ to investigate more exotic superfluid phases. 
On the theoretical side vortex core states in superfluid 
fermion mixtures have been extensively studied for mass 
and population balanced mixtures [ll|, [l^, [H, [13, flU . 
In these works it was predicted that the local density of 
fermions is slightly depleted at the vortex core for weak 
interactions, and that the density depletion increases 
with increasing interaction strength. Furthermore vortex 
core states of mass balanced but population imbalanced 
fermion mixtures have recently been analyzed [TgI . [iTj . 
showing that the unpaired (excess) fermions occupy the 
core states. However these predictions have not been 
experimentally observed since a probing technique anal- 
ogous to the scanning tunneling microscopy is still lack- 
ing. 

Arguably one of the current frontiers of ultracold 
atom research is the investigation of two-species mass 
imbalanced fermion mixtures (e.g. ^Li and '^"K, ^Li 
and ®^Sr, or '*°K and ^^Sr) with or without a pop- 
ulation imbalance, due to their greater potential for 
finding exotic phases [H, [l^. Such mixtures are cur- 
rently of interest to many communities ranging from 
atomic and molecular to condensed and nuclear matter 
physics, and recent analysis of the ground state phase di- 
agram have shown quantum and topological phase tran- 
sitions 0, iH, m, [2I, m m, nil. in this manuscript, 
motivated by the very recent experimens on ^Li and ^''K 
mixtures [isl . [Toj . we analyze the vortex core states of 



superfluid Fermi-Fermi mixtures with unequal masses as 
a function of two-body binding energy. In particular we 
use the Bogoliubov-de Gennes (BdG) formalism to study 
a population balanced mixture of ^Li and ^'^K atoms at 
zero temperature. We find that the vortex core is mostly 
filled with the light-mass (^Li) fermions and that the core 
density of the heavy-mass (^"K) fermions is highly de- 
pleted. This is in contrast with the one-species (mass 
balanced) mixtures with balanced populations where an 
equal amount of density depletion is found at the vortex 
core for both pseudospin components [ll|, [H, [H, [13, [l^ . 

We achieve these results by using the following 
Hamiltonian density H{x) = X^o- ^ct(-^)^<^('')^<^(''') ~ 
g"^^ {x)^{x), which describes two-component fermion 
mixtures with attractive {g > 0) and short-range in- 
teractions. Here h = ks — I, and ipaix) 
are the Grassmann field operators corresponding to 
creation and annihilation of pseudospin a fermions at 
position r and time t {x = {r, r}). Furthermore 
K,{r) - -VV(2M,) - fi, and *(a;) = 4,^{x)i:^{x), 
where Ma- is the mass and /io- is the chemical poten- 
tial of (T fermions. In the mean-field approximation for 
the superfluid phase the resultant Hamiltonian can be 
diagonalized via the Bogoliubov-Valatin transformation 
V'^tW = E„K,<T(r)7„,a-s<T<,^(r)7l where u„,^(r) 
and fn,o-(r) are the amplitudes and 7^ ^ and 7„^cr are the 
operators corresponding to the creation and the annihi- 
lation of pseudospin a quasiparticles, and S| = -1-1 and 
S| = —1. Then the single particle Green's function ma- 
trix (in Nambu pseudospin space) can be written as 



G{r,r' ■,iLUe) = ^ 



(1) 



where u)e = {2i + 1)ttT is the fermionic Matsubara fre- 
quency, £ is an integer number and T is the temperature. 
Here V'n.o-Cr) and en,a- > are the eigenfunctions and the 
eigenvalues of the BdG equations 



A*(r) 



A(r) 
-Kliv) 



Vn.air) = Sa£n,aVn.a{r), (2) 



where (pn,a{r) is given by (pn,^Hr) = [m* -|-(r), j^(r)] for 
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the t and (pn,i''{r) = [w„,f (r), -M„_|(r)] for the i eigen- 
values. Since the BdG equations are invariant under the 
transformation Wnj(r) — > u* |(r), ■u„_|(r) — ^(r) 
and en. I — > — Cri.ti it is sufficient to solve only for 
M„(r) = Un^^ir), Vnir) = w„,x(r) and e„ ee e„,| as long 
as we keep all of the solutions with positive and negative 
eigenvalues. 

In Eq. ^ ^(r) is the local superfluid order 
parameter defined by A(r) — (7('0f (r)-!/)j^ (r)) — 
— gT^^ G|_|(r, r; iw^), which after evaluat- 
ing the frequency sum leads to A(r) — 
-9J2n,aS'^UnA^)K-A^)fisaen,a)- Here (...) is a 
thermal average and f{x) — I/[exp(a:/T) + 1] is 
the Fermi function. Using the symmetry of the 
BdG equations this equation can be written as 
A(r) = -5S„"n(r)^^«(r)/(en)- We also relate the 
interaction strength g to the two-body binding en- 
ergy eh < of an I and a J, fermion via the relation 
1/g = (IM) Ek l/(<^k,T + ek.i - Cfc), where A is the 
area of the sample and et.cr — fc^/(2Mcr) is the kinetic 
energy. This leads to g — 47r/[Mrln(l — 2ec/eb)], where 
Mr — 2M|M|/(Mf -f- All) i^ twice the reduced mass of 
an "f and a J, fermion and Ec is the energy cutoff used in 
the k-space integration. The order parameter equation 
has to be solved self-consistently with the number equa- 
tions Ncr = J drna{r) to determine fi^, where ^^(r) = 

is the local density of a fermions. After eval- 
uating the frequency sum this relation leads to 
na{r) = E„[|wn,<T(r)P/(e«,<T) + \vn,a{r)\^ f{-en^-a)], 
which can be written as n|(r) = J2n ^^'^ 
= J2n l^n(i')P/(-e„) by using the symmetry of the 
BdG equations. Having discussed the BdG formalism, 
next we analyze the self-consistency equations for a 
single vortex. 

In particular we consider a two-dimensional homoge- 
nous disk of radius R such that the local order param- 
eter can be written as A(r) = A(r) exp(— ikS), where 
r — {r, 0) are the polar coordinates and k is the vor- 
tex winding number. This choice is due to numeri- 
cal reasons and we do not expect any qualitative dif- 
ference between our results and the three-dimensional 
ones. Then the normalized wave functions are of 
the form u„(r) = exp(im6')/-\/27r and ti„(r) = 

Vn,m{j') exp[z(TO-|-K)0]/-\/27r such that the BdG equations 
can be solved separately in each subspace of fixed angu- 
lar momentum m [ll[. We further project the radial 
wave functions u„,m(r) = Y,j Cn,j<l>j,m,{r) and w„,m(r) = 
dn,j4>j.m+K{'r) onto a set of orthonormal Bessel func- 
tions 0j,m(r) = V2Jm{aj,mr / R) /[RJm+i{aj,m)], where 
the argument aj,m is the jth zero of Jm{x). This pro- 
cedure reduces BdG equations given in Eq. ^ to a 
2jmax X 2jmax matrix eigenvalue problem [ll| 



E 



I ,m 



AO -3 



= ^«E 



,(3) 



„/(2Mo.i?^) — Hcr\5j.j' is the diagonal and A^ 
JrdrA(r)(l>j^„i{T)(f>j',m+K.{'r) is the off-diagonal element 
where Si^ is the Kronecker delta. Furthermore the order 
parameter equation reduces to 



2tt 



and the local density equations reduce to 



rM = E -^^^<t>jMr)<Pf,rn{r)fien), (5) 



iw = E 



^,~Ar)^,,~{r)f{-e^), (6) 



2tt 3:1^^ '^]',m 



where m = m + k. Notice that the total numbers of t 
and i fermions are given by N-^ = J2n,m,j '^l.jfi^n) and 
= Z^n.mj dl .jf{-en), respectively. We emphasize 
that these mean-field equations can be used for all val- 
ues of g but they only provide a qualitative description 
of BCS-BEC crossover at zero temperature (T = 0) as 
discussed next. 
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FIG. 1: (Color online) Local density of a fermions ncr{r) [in 
units of fc|./(27r)] versus radius r (in units of l/kp) is shown 
for a population balanced mixture of ^Li and atoms. Here 
\eb\ — O.lep (solid lines) and \eb\ = 0.2eF (dotted lines). The 
inset shows the local order parameter A(r) (in units of ef) 
versus r for the same parameters. 



In this manuscript we discuss a population balanced 
mixture of ®Li and ^"K atoms, where N-^ — and 
Af| = 0.15Af|. For this purpose we set a large energy cut- 
off ec = Oep, and numerically solve the self-consistency 
Eqs. ([3]), (HI), (O, and © for a singly quantized vor- 
tex with K = I at T = 0. Here ep — kp/{2AIr) is a 
characteristic energy scale where kp is the Fermi mo- 
mentum corresponding to the total density of fermions. 
We also choose R — 25 /kp as the radius of the two- 
dimensional disc, and jmax = 25 and |TO|maa; = 50 as 
the maximum quantum numbers. Notice that the bulk 
density = N^/A is given by ni = ni = kp/{ATr) for a 
population balanced mixture. As one may expect pres- 
ence of a single vortex can not significantly effect the bulk 
parameters. Therefore, to simplify the numerical calcu- 
lations, we first solve /Xo- and |Ao| self-consistently for a 
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vortex-free system, and then use these solutions as an 
input for our vortex calculation. Here |Ao| corresponds 
to the bulk value of A(r). Since our vortex calculation 
is not fully self-consistent, and Ni turn out to be 
very close |7V| - 10~^{N^ + N^) but not exactly 

the same. Furthermore, for mass balanced mixtures, we 
checked that this procedure gives results that are in qual- 
itative agreement with the earlier works on population 
balanced [ll|, [T^ 13 . [Tsj as well as population imbal- 
anced svstems~jl6l. Il7l|. 

In Fig. [T] we show n^{r) for a population balanced mix- 
ture of ^Li and ""^K atoms. We also show A(r) as an in- 
set for the same parameters. When g is small such that 
|ef,| ^ ep, the vortex core is mostly filled with the light- 
mass (®Li) fcrmions while core density of the heavy-mass 
(''"K) fermions is highly depleted. This is because the 
bound state energy spectrum is discrete with a small but 
nonzero separation ^ |Aop/eF (discussed below), and 
thus only the light-mass fermions can occupy these states 
at the core since fj,Li > Mis' due to the mass difference. 
While this is in sharp contrast with mass and population 
balanced mixtures where the bound states are unoccu- 
pied and equal amount of depletion occurs for both a 
fermions [HI, [13, [13, Hi EBl , it is similar to mass balanced 
but population imbalanced mixtures where the core is 
filled with excess fermions due to their higher chemical 
potential [l^, [13 • However, as shown in Fig. [1] local 
density of the light-mass fermions as well as that of the 
heavy- mass one deplete more with increasing |e;,|, which 
is qualitatively similar to that of mass balanced mixtures 
where the local density depletion also increases with in- 
creasing g [13, [3, [iB| ■ This is because the energy separa- 
tion between the bound states increases with increasing 
|Ao| which makes them less occupied. To further under- 
stand these peculiar density depletions, next we analyze 
the single particle density of states for a fcrmions as well 
as the spectrum of energy eigenvalues. 

At r = the local single particle density of 
states for a fermions is defined by Dcr{r,uj) — 
— (l/7r)Im[limg^Q+ G'cr^CT(r, r; ^^^f clJ + ie)]. This leads 
to D^{r,uj) = T,n{Wn,air)\'^S{uj - en.a) + \vn,a{r:)\^5{uj + 
en,-iy)\ where 6{x) is the delta function, and it can 
be written as D^{v,lu) = X^n ~ ^n) ^^"^ 

Dx(r,w) = Y.n |i'n(r)P'5(w -I- e„) by using the symme- 
try of the BdG equations. Then the overall single parti- 
cle density of states is found by Dcr{i-o) = J drD^{v,uo), 
which for a single vortex reduces to D^{lo) = 
En,mj - e«) and Di{uj) = E„,mj dl^j^iuj + e„). 

We use a small spectral broadening (O.Olep) to regularize 
these delta functions. 

In Fig. [2] we show D„{uj) for a population balanced 
mixture of ^Li and ''"K atoms. We also show the spec- 
trum of energy eigenvalues e„ as an inset for the same 
parameters. Similar to mass balanced mixtures [ill , [l^ . 
[isl . [3, [l3|, the positive and negative energy spectra 
are connected by a single branch of discrete Andreev- 
like bound states. The visible discreteness of the con- 
tinuum spectrum shown in Fig. [2] is a finize size ef- 




FIG. 2: (Color online) Single particle density of states for a 
fermions Dcr{uj) (in units of I/ef) versus energy uj (in units 
of ef) is shown for a population balanced mixture of ®Li and 
^''K atoms. Here |e(,| — O.leF. The inset shows the spec- 
trum of energy eigenvalues e„ (in units of ef) versus angular 
momentum m for the same parameters. 



feet and the spectrum becomes continuous only in the 
thermodynamic limit {kpR — > oo), while the discrete- 
ness of the bound states is insensitive to the system size 
since these states are strongly localized around the vor- 
tex core. However the energy spectrum is asymmetric 
around e„ = with less bound and continuum states 
for e„ > 0, which is due to the broken pseudospin sym- 
metry since the masses {mx > Tnn) and therefore the 
chemical potentials {^iu > ^J■K) are different for t and J, 
fermions. This is in sharp contrast with mass and pop- 
ulation balanced mixtures where the energy spectrum is 
symmetric [3, [l3| . 

Furthermore the energy spectrum is qualitatively dif- 
ferent on the positive and the negative to sides, and it is 
very illustrative to make an analogy between the energy 
spectrum of mass imbalanced mixtures shown in Fig. [2] 
and that of the mass balanced mixtures [l3,[iB| • For mass 
and population balanced mixtures energy spectra that 
are qualitatively similar to the negative (positive) side 
with many (few) bound and continuum states occur for 
small (large) values of |eb|, leading to low (high) density 
depletions at the vortex core. This analogy suggests that 
the local vortex core density of the heavy-mass fermions 
should deplete more than that of the light-mass fermions 
since the density of states for heavy-mass fermions is 
higher (lower) for positive (negative) ui. We also find that 
the bound state contribution to n^r (r) is a nonmonotonic 
function of r with a maximum at an intermediate dis- 
tance r = r*. This nonmonotonic contribution is due to 
the strongly localized quasiparticle amplitudes that are 
associated with the bound states, which also give rise to 
Friedel-like A(r) oscillations around the vortex core in 
the strict BCS hmit [uj, [13, [H, [ll- Therefore the un- 
equal density depletions shown in Fig. [T]are purely den- 
sity of states effects arising from the asymmetric energy 
spectrum shown in Fig. [21 Next we analyze the density 
as well as the velocity of the superfluid fermions. 

The quantum mechanical probability cur- 
rent operator for a fermions is given by 
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J<,(r) = [l/(2M,i)][i^t(r)VV'a(r) - H.c] where 
H.c. is the Hermitian conjugate. Therefore the 
focal current density Jcr{^) = (Jcr(r)) circulat- 
ing around a single vortex becomes J| (r) = 
[l/(2Mt»)]E„[<(r)Vu„(r)/(e„) - H.c] for the t 
and Ji(r) = [l/(2Mx*)] E„[««(r)V<(r)/(-e„) - H.c] 
for the I fermions where we used the symmetry of 
the BdG equations. These relations can be written 
as Jo-(r) = no(r)vo-(r)/2, where no(r) is the local 
superfluid density and Vo-(r) — K9/{2Ma-r) is the 
local superfluid velocity. Therefore J(T(r) is along 
the 9 direction, and for a single vortex it is given by 
JtW = [l/(27rMTr)]E„.,„mE,c„,,0,,™(r)]2/(en) 
for the r and Jx(r) = -[l/(2^Mxr)] + 
dn,]4>],m+K{r)]'^ f{-en) for the I fermions such 
that no{r) = {4A'U/ K.)rJ„{r). 
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FIG. 3: (Color online) Local current density J|(r) [in units 
of Mi/k%\ versus radius r (in units of l/kp) is shown for a 
population balanced mixture of ^Li and atoms. Here 
\et\ — O.lep (dashed line), |ei,| — O.Sep (solid line) and 
\eb\ = O.Sep (dotted line). The inset shows the local super- 
fluid fermion density no(r) [in units of fc|-/(27r)] versus r for 
the same parameters. 

In Fig. [3] we show J|(r) for a population balanced mix- 



ture of ^Li and ""^K atoms. We also show no(r) as an 
inset for the same parameters. The bound states have 
positive (paramagnetic) and the continuum states have 
negative (diamagnetic) contribution to J| (r) . This leads 
to a nonmonotonic J| (r) which first increases as oc r and 
then decreases as oc 1 /r. The latter behavior is due to the 
saturation of no (r) for long distances away from the vor- 
tex core. Therefore a maximum peak current occurs for 
all values of |ef,| at some distance Vc away from the vor- 
tex core. However the value of this peak current increases 
until jefjj ~ 0.3ep and then decreases for higher values of 
|ef,|. Since a two-body bound state exists even for an 
arbitrarilly small 5 > 0^ in two dimensions, we empha- 
size that this nonmonotonic evolution is not due to the 
occurrence of a two-body bound state threshold (diver- 
gence of the two-body scattering length), as previously 
suggested for mass and population balanced mixtures in 
three dimensions ■ We believe that it is related to the 
nonmonotonic evolution of the coherence length which 
can be easily extracted from no{r). This is qualitatively 
consistent with the recent experiments involving mass 
and population balanced mixture of ®Li atoms, where a 
pronounced peak of critical velocity has been observed 
on the molecular side of the stron gly interacting regime 
in one-dimensional optical lattices [27[. 

In summary we analyzed the vortex core states of pop- 
ulation balanced ^Li and mixture at T = as a 
function of two-body binding energy. We found that the 
vortex core is mostly occupied by the light-mass (^Li) 
fermions and that the core density of the heavy-mass 
(""^K) fermions is highly depleted. This is in contrast 
with mass and population balanced mixtures where an 
equal amount of density depletion is found at the vortex 
core for both pseudospin components. 

This work was partially completed in Georgia Tech, 
and the author gratefully acknowledges discussions with 
C. A. R. Sa de Melo and P. S. Julienne. 
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